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ABSTRACT 


A stochastic analysis of the Supply Corps officer 
manpower system is presented. The analysis is based on the 
concept of diffusion approximations. The state of a rank 
is represented as the superposition of a deterministic process, 
and a random noise variable, both of which are derived 
utilizing a set of formulated recursion equations.  Numeri- 
cal results are presented for a 3-rank system and compared 
against a simulation of a Markovian manpower model. 

A Diffusion Approximation analysis is performed utilizing 
United States Navy Supply Corps data, in an effort to 
determine the probability of maintaining each rank, any 


two ranks, and/or the entire system between specified limits. 
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I. INTRODUCTION 


This paper attempts to examine manpower problems uti- 
lizing stochastic techniques. The point of departure is the 
deterministic theory of control given in Bartholomew [Ref. 1], 
and the object is to see how concepts and strategies origi- 
nating in that context perform in a stochastic environment 
of the kind likely to be encountered in practice. Conclusions 
reached, both in this thesis and in Bartholomew [Ref. 3], 
suggest that, although useful, deterministic analysis is 
incomplete since it fails, for example, to provide required 
information concerning the variance of individual ranks. 
Without such information, manpower planners may attempt to 
"fine tune" the system by increasing or decreasing promotion 
or recruitment rates when in reality the fluctuations are 
due to the uncertainty in the system. 

This study concentrates its efforts on one particular 
manpower system--that being the Supply Corps of the United 
States Navy. An effort is made to develop a stochastic model 
that provides the capability for simple computation. In 
addition, it will give the mean number of individuals in 
each rank, as do deterministic models; and as added features, 
will provide second order moments and distribution of force 
information. Armed with these key elements we can then 
attempt to answer questions related to the maintenance of 
desired force structures and with the ability to attain such 


structures. 





The numbers of people in particular ranks change over 
time as a result of attrition, promotion and recruitment. 
Some of these flows are subject to management control while 
others vary due to a multitude of circumstances and may best 
be modelled stochastically. Factors such as the need to offer 
adequate career prospects, while at the same time, maintain- 
ing prescribed force levels are the manpower planners chief 
goals. A problem arises in that the planner must determine 
whether these goals can be achieved simultaneously, and, if 
so, how best to achieve them. 

The manpower control problem has received a good deal of 
attention in recent years--see, for example, Bartholomew 
[Refs. 1,2,3], Davies [Refs. 4,5] and Grinold and Stanford 
(Ref. 6]. 

Although most of the work on manpower control has been 
based on stochastic models, the analysis has been determinis- 
tic in the sense that only expected values for attrition and 
promotion rates are used. More realistically, attrition and 
other flows do vary randomly and are therefore better repre- 
sented by random variables. Recent work which does attempt 
to incorporate randomness into manpower control are Bartholomew 
[Refs. 2,3] and Kim [Ref. 7]. 

In this thesis we will assume that only promotion rates 
and recruitment rates are subject to control. While it may 
be true that attrition can be dramatically increased, uti- 


lizing techniques such as "Reductions in Force", it is the 





problem of being unable to hold down attrition rates that 
plagues the present day Navy. As a result, the model developed 
in this thesis utilizes, as an input parameter, the attrition 
rate and is capable of determining, in a deterministic sense, 
the promotion and recruitment rates required to maintain the 
system at some specified level.  Henceforth, the concept of 
being able to maintain a specified force structure will be 
referred to as maintainability; i.e., a structure is main- 
tainable if and only if there are promotion and recruitment 
rates such that for each rank in the system, the expected 
number of individuals who depart during each time period is 
equal to the expected number of individuals who arrive during 
that same time period. The majority of this thesis will deal 
with the maintainability of desired force structures, but it 
turns out that this topic cannot be treated adequately with- 
out considering the question of whether those structures can 
be attained from present levels. 

Since for all practical purposes, recruitment into the 
Supply Corps takes place only into the lowest rank, promotion 
control alone is left as a means of controlling the force 
distribution. Control by promotion has the disadvantage that 
fluctuations in the promotion rate are liable to create 
internal stresses within the organization. On the other 
hand, it was shown in Bartholomew [Ref. 1, Chapter 4], that 
in a deterministic environment, promotion control is capable 
of maintaining a much wider range of structures than recruit- 


ment control alone. 





Having discussed the various areas of consideration, the 
general scope of this paper will be to present the following: 

A. A basic deterministic Markovian model will be 
derived in Chapter II, to be used as a benchmark with which 
to compare the stochastic version. 

B. Realizing that a deterministic representation is not 
representative of the random factors that affect most pro- 
cesses, a more realistic stochastic model will be derived 
in Chapters III and IV utilizing a Diffusion Approximation. 

C. Results of a simulation of a stochastic Markovian 
manpower model will be presented in Chapter V and compared 
to results of the Diffusion Approximation as developed in 
this thesis.  Shortfalls of the approximation will be dis- 
cussed and advantages of the approximation pointed out. 

D. A comparison of the deterministic and stochastic 
models will be made in Chapter VI. In particular, the corre- 
lation between adjacent ranks, and its effect on the ability 
to maintain desired force structures will be discussed. 

E. Areas of possible further study will be briefly 
mentioned in Chapter VI. An effort is made to identify 
possible research topics, whose results, when used in conjunc- 
tion with the model derived in this thesis, may reflect to 
an even greater degree the stochastic nature of a manpower 


system. 





II. DETERMINISTIC ANALYSIS 


A. A DETERMINISTIC MANPOWER MODEL FOR THE SUPPLY CORPS 

Before one attempts to model any system, one must be 
familiar with its everyday operations. In an effort to 
compare the results of the approximation against a benchmark, 
we will use our familiarity of the system to make several 
educated assumptions as to the values of promotion and attri- 
tion rates within the supply Corps. A compilation of these 
educated guesses, will define a Markovian transition matrix; 
the elements of which are the fractions of individuals in 
each rank that stay in the rank or are promoted in each time 
period. This matrix will then be compared against one where 
historical data is used as a basis for determining the ele- 
ments of the Markovian transition matrix. Comparison between 
the two will answer the question of whether organizational 
assumptions made are reasonable. 

The organizational assumptions that we will make are as 
follows: 

1. We will assume that the fraction of Supply Corps 
officers with a particular number of years in a specified 
rank remains almost constant; that is, we will assume that 
ene fraction of the Supply Corps force that are second year 
Lieutenants or third year Lieutenant Commanders, etc., 
remains the same. 

2. Under present promotion patterns, promotion from 


Ensign (0-1) to Lieutenant Junior Grade (0-2) occurs at the 


10 





2 year point in an officer's career. In addition, promotion, 
for all practical purposes, is assured for all Ensigns com- 
pleting their second year. It is therefore assumed in light 
of assumption l, that approximately 50% of all O-l's are 
completing their second year and will therefore be promoted 
mo, O-2. 

3. Similarly, it takes two years to transition between 
ranks O-2 and Lieutenant Commander (0-3); however, only 95% 
of those in the second year in rank O-2 get promoted. There- 
fore, this implies that 50% of all O-2's are in the promotion 
zone, and 95% of those eligible are promoted. This informa- 
tion yields the promotion fraction equal to 0.475 (0.50 x 0.95). 


4. Using similar logic the following table values are 


computed: 
Years Required Assumed per- Assumed Promotion 
in grade prior centage of promotion fraction of 
to being in total grade rate for all members 
zone in zone in zone in grade 
Ensign (O-1) 2 0.5000 1.0000 0.5000 
Lieutenant 
Junior Grade 
(0-2) 2 0.5000 0.9500 0.4750 
Lieutenant 
(0-3) 5 0.2000 0.9000 0.1800 
Lieutenant 
Commander 
(0-4) 7 0.1429 0. 8000 0.1143 
Commander 
(0-5) 5 0.2000 0.6000 0.1200 


LI 





Τη order to complete the transition matrix, it is necessary 
to determine reasonable attrition rates.  Attrition rate data 
was obtained from historical data files at the Naval Military 
Personnel Command.  Attrition rates cannot be forecasted with 
great accuracy, due in great measure to their correlation to 
many economic and social factors. Thus it appears reasonable 


to use a simple moving average as a predictor. We will use 


Vi (t-1) t Vi (t) 
2 


as a forecasted value for the current attrition rate for 
Stage 1 at time t+l. 

This average Bora to be a better forecasting technique 
than that of simply using the last attrition rate. A reason 
for not considering more terms in the attrition forecast, 
is the fact that we want the forecasted attrition rate to 
reflect only the recent past and not to consider forces which 
may have affected the manpower system at one time but are no 
longer present. 

The remaining element required in the transition matrix 
is the continuation rate. Since the sum of the continuation 
rate, the attrition rate and the promotion rate must equal 
one, it is a simple matter to compute the remaining elements 
of the transition matrix. It should be mentioned once again, 
that these assumptions are not related to any historical 
data, except for attrition rates and existing "in zone" 
promotion rates. The resulting matrix is computed merely to 


compare against the matrix derived in the next subsection. 
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B. PROMOTION RATES REQUIRED TO MAINTAIN THE SUPPLY CORPS 
AUTHORIZED FORCE SIZE 


The objective of this subsection, is to find out if the 
force levels which the supply Corps has been authorized, 
are in fact maintainable. The authorized force levels men- 
tioned are included in Appendix A, and are obtained from the 
Officer Programmed Authorizations, Fiscal Years 1979-1983 
(Ref. 9]. In an effort to work with a single force level, 
the average force level was computed for the five years men- 
tioned above; and it is this level that we attempt to maintain. 
It should be noted that over the 5 years under consideration, 
neither the individual ranks nor total force strength vary 
significantly, so the average force appears to be a reasonable 
force level to maintain. 

If we determine that the desired force can be maintained, 
we will then have a transition matrix which can be compared 
against the matrix Pi derived in subsection II.A. This will 
give us an idea of how drastically the system would have to 
change in order to maintain the desired force levels. 

Before we proceed, a few definitions for the above: 


Q = a six component row vector whose elements 
are the number in each rank. 


P z a 6x6 matrix whose elements contain promo- 
tion and continuation rates. (See matrix 
P,, subsection II.A for the general form 

ot this matrix.) 


V = a column vector whose elements contain the 
; appropriate attrition rates for each rank. 
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Ros Ἶποικεαιααιαποης vector, which directs the 
~ flow of recruits into the proper ranks 
(fraction of recruits for each rank). 

In the Supply Corps, the recruitment vector is not gen- 
erally amenable to control, i.e., recruitment for the most 
part occurs into the lowest rank, i.e., R= [1,0,0,0,0,0]. 

In such cases we wish to know what degree of control can be 
exercised by the promotion rates. For a fixed size organiza- 
tion, what we wish to do is to find a transition matrix P, 
which allows us to maintain our desired structure Q. This 


implies that P must satisfy 
Omran or + 0V R 18:1) 


where Qvi R yields the number of recruits necessary to 
exactly equal the number of attritions at each time period. 
As was noted in subsection II.A, the matrix P is super- 
diagonal in form, i.e., non-zero elements only in the main 
diagonal and the diagonal above it. Bartholomew [Ref. 1] 
shows that when P has the super-diagonal form, equation (II.B.1) 
admits a unique solution which may be tested to see whether its 
elements are non-negative. If the matrix P is found to be 
non-negative, then the structure Q can be maintained. 

As was mentioned previously, while we assume that 
recruitment takes place only into the lowest rank, we recog- 
nize the fact some recruitment does occur in ranks O-2 and 


O-3. However, the number recruited into these ranks is small 


15 





and hence will be neglected. Also, the attrition rates 


are as considered previously, i.e., 


a2 727027 355,0.1121,0.0390,0.0790,0.1622]. 


We note that equation (II.B.1) is equivalent to 


and, at this point it becomes more convenient to abandon 
the matrix notation temporarily and to proceed with the com- 


ponents as does Bartholomew [Ref. 1]: 


k 
τ; l Qiy 


al NINE 
1=1 
k 
E PE 
k 
Q4P56 t Q4Pgg = Ὡς 7 ἃς ;28iVi 


(The above are components of the respective vectors or 

matrices mentioned at the beginning of this subsection.) 
Eliminating Pii using Pii = 1 - Vi Ξ ΤΙ Gin ον 

k-1) and using our special structure for the vector R yields 


the following equation: 


l6 





K 
Po; = ο 70. απ ο) 
τι jeitl RE RE 


In other words, equation (II.B.2) says that the proportions 
required to be promoted from grade i must equal the number 
leaving from all grades above grade i divided by the size 
of grade i. 

For example, given a point estimate for the attrition 
rate, and given that the promotion rates for O-6's is essen- 
tially zero; it is possible to compute the number of O-6's 
who depart the system. Given this value, one can then 
compute the number of O-5's who must be promoted each year 
to exactly maintain the required level of O-6's. The process 
continues in this fashion until we compute the required num- 
ber of new accessions required to maintain the O-1 Rank. It 
is in this way that promotion rates can be calculated to 
exactly maintain the desired structure (assuming that we are 
at the structure we wish to maintain). If we are at some 
other force level and we wish to attain our goal, similar 
calculations can be done, but review must be done to ensure 
that the resulting promotion rates are within acceptable 
bounds as specified by organization policy. 

For example, suppose we wish to maintain the following 


structure: 
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Total Force Level 5 3955 





Rank Strength Attrition Rates 
dl 621 0.0272 
2 625 021355 
2953 995 0.1121 
O-4 863 020390 
05 617 0.0790 
O-6 234 ο το” 


We first compute the number of O-6's lost. This is equal to 
(234) (0.1622) = 37.9548. To maintain the desired structure 
37.9548 O-5's on the average must be promoted. Therefore, 
the promotion rate required is (37.9548) “ (617) = 0.0615. 
An important note to make, is that promotion rate is used 
in relation to the total rank strength, to determine what 
the rate would be in regards to "In the promotion zone" 
strength, one must first determine the number to be promoted, 
and then divide that number by the number "in the zone". In 
our example, promotion rate related to "In the zone" would be 
37.9548 — (number of O-5's in zone). The other promotion 
rates can be similarly computed, and as a final result we 
obtain the matrix P5. 

Comparing this matrix with P1 reveals that the reasonable 
guesses for the existing system are very close, in most cases, 
to the rates necessary for maintaining the desired system on 


the average. This implies that very little has to be done in 
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way of changing promotion rates if we simply want to get 
to and maintain the desired levels. It should be stressed, 
that once a specified force level has been shown to be 
maintainable, it can be attained from any starting force 
level as was stated in Ref. 1. In the Supply Corps, present 
rank levels are not equal to desired rank levels. However, 
if promotion rates were modified to those suggested by Ρο» 
we would eventually reach a system of manpower stocks whose 
expected value in each rank would equal the desired levels. 
A technique that suggests a method to speed up the process 
of convergence from present levels to desired levels will be 


discussed in a later chapter. 
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ΤΙ απ ΡΕΡΕΕΟΡΤΟΝ MODEL 


The system we wish to model is as depicted below: 


— 





Ὡς (t) 


9 ο - Qc (€) 


Ve (t) 


The variables of interest are as follows: 


a(t) 


u; (t) 
v; (t) 


ο, (5) 


the force available for recruitment, 
assumed to be large. 


the accession rate at time t 


the promotion rate for all members in 
rank i at time t 


the attrition rate for all members in 
rank i at time t 


the Random Variable which is the force 
level in rank i at time t. 


Other variables of interest to be incorporated into the 


model are: 


q, (t) 


the expected fraction of the total force 
in rank i at time t (computed as fraction 
of original total force level) 


E[Q; (t)] 


or 
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Wi (t) = Promotion random noise variable, 
distributed N(0,1) (when i = 0, 
Ww P(t) relates to accession) 

wW,” (t) = Loss random noise variable, distributed 
Π.Ο, 

X; (t) = standardized random noise variable. 


The probability that a member of rank i gets promoted 
at time t is equal to ui: (t); this implies that the proba- 
bility that the member does not get promoted is equal to 
l-u;(t). We will assume that whether or not an individual 
is promoted is independent of others in the system at time 
t. Thus, the number promoted will have a binomial distribution 
with expected value u, (t)Q, (t) and a variance of 
u, (t) [L-u, (t)1Q, (t). 

Using similar reasoning the number of attritions in rank 
i at time t is assumed to have a binomial distribution with 
an expected value of vi (t)Q, (t) and a variance of 
v; (t) I1 7 vi (£)1Q, (6). 

In rank 1 we have a special case, since we must consider 
the accession rate as a function of the total force level. 
We will assume that there is a large population of potential 
recruits of size Qu We assume the number of people that 
enter the system at time t has a binomial distribution with 
mean a(t)Q. and a variance of a(t) [1 - a(t)]Q0,- In what 
follows we will usually take Qo to be the initial force 


strength. 
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If one deals with a system which is sufficiently large, 
one can invoke the Central Limit Theorem. Therefore if Q 
is "large", this would imply that Q; (t) is "large" and 
therefore the distribution of number promoted can be approxi- 
mated by a Normal Distribution with mean u; (t)Qi (t) and a 
variance of αι (5) [1 - u, (t)]0, (t).. Similar results apply to 
the number of attritions in each rank. 

Now, recall if X is distributed normal with mean u and 


variance ix - πο. Then 
aa). 


Let 


then 


E 


Qi; (t)ui; (t) + "Qi (t) u, (€) (1 - u, (t)] πι (t) 


3 N(Q; (t)u; (£),Qi (t)u; (€) [12 7 ui (€) T) 
which is a normal distribution with mean and variance equal 


to the mean and variance specified for the binomial distribu- 


tion of the number promoted at time t. 
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Similarly, 


| - b 
Q; (t) v; (t) + VQ, (t)v, (t) [1 v,(t)] W, (t) 


is a normally distributed random variable with mean and 
variance equal to the mean and variance specified for the 
binomial distribution of the number attrited at time t. 
For each interval of time (t,t+1] we can define the 
change in rank i, i.e., Q; (t+1) - Q; (t), as a combination 


of normally distributed random variables. 


a a o τε P 
DE O te) Ww..." (t) 


H i-l 


Eiei- 


(promotions into rank Q;) 


A EE reacio a ut 
u; (t)0, (t) u; (t) [l u; (€) JOi (t) Wi (t) 


(promotions out of rank ο.) 


- v, (t)Q,(t) - w, ilv, la, E W, (5) 


(attritions out of rank Q;) 


From this relationship, Q; (t*1) can be solved for in terms 


of Qi (t). Hence, we assume for i» l, 
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= πμ... Πο. (ey 5 
Q; (t+1) u; mo, (8) + Yu, (8) [1-u,_, (t)]0,_,(t) W,_, (€) 


- ¿(Eu (TO (e W; (t) 
- „rw, TE) II -v,te)To,te) w,(t) 


+ [1 - u, (t) - v, (t)] Q; (t) (oT 


where (W, (t), i = 0,1,...,r-1) and (W,"(t), i - l,...,r) 
are independent random variables, having a normal distribution 
with mean zero and variance equal to one. In addition all 
the W random variables are independent of (Qi (t), LS EIS 
inet 

ο, (5) - 0,9, (t) 


X; (t) = SSS σετ 2) 


Lo 


It Qo is sufficiently large, then X(t) = [X (EI, X, (EI) 1... ,X, (8) ] 
should have approximately a multivariate normal distribution 


with mean zero. 
From equation (III.2) we have 
Qi (t) = Qog; (t) + “O, X; (t) jeher ZT RR) 
In Appendix B we replace the variables Q; (t) in equation 


(III.1) by the expression in (III.3) and then take the limit 


as 0, > «e. This procedure yields recursive equations for 
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(q; (t),1 ana (X. (t),i zum εἰ The 


l,...,r)] are constants and (X; (t),i —p Ir 


(q; (t), i 
is a multivariate normal random variable with mean zero. 
Hence in light of (III.3), Q(t) will also have a multivariate 
normal distribution, at least in the limit as o > o. In 

the next chapter we will present a recursive scheme to 


evaluate q; (t), Var [X, (t) ] and ο νο ου 
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IE CURSTEVE SCHEME 


Only results will be presented in this section.  Deriva- 
tions of the results will be presented in Appendix B of this 
paper. 

The following are definitions required for the recursions 


which follow: 


Var [X, (t) ] m, (t) 


TT 
3 
ct 


E A 


RESULT l: A RECURSIVE SCHEME FOR COMPUTING αι (9), E T 


i = 1: q, (t*1) a(t) -* [1 - u, (t) - vi (t)] q; (t) 


mE I: q; (t*1) u¡-7 (619,1 (t) + [1 - u, (t) - v, (t)] q; (t) 


RESULT 2: A RECURSIVE SCHEME FOR COMPUTING Var [X,(t)], to eens 


isl: m(ttl) = [L-u,(t)-v,(t)]* m(t) 
SS E) 


+ tu, (t) [1 - u, (t) ] + v, (t) [1 - v, (t)15q, (€) 
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i> dl: m (ttl) = fu, Q(0)]^ m, ,(€) 
2 

* [2-ui(t) "v; (€)] mj (t) * ui (1-u,_,(t)lg,_,(t) 

+ (Qu, (E) [7 u; (€) ] e v, (E) [1 7 v, (€) 1) q, (6) 


+ 2u; j(t) [1 - u, (t) - vi(t) lmi 1, i © 


RESULT 3: A RECURSIVE SCHEME FOR COMPUTING COV [X, (EI ^ (eus 
e=- l2; 


W l and j = 1+1: 


mya Ξ α, γ(Έ)α; (Ε}πι. 4 (8) 


+ uj (E) [1 o5 (6) = v(t) Im, (6) + uj (t) 


+ [1 - u; (t) - vi (t)] LEES ANN 


- u; (t) [1 - u; (t) Iq, (t) 


EE l and j #7 i+l 


Mas E = A O q (Em, _ los -1 (6) 
+ u; (t) ΠΡ u. (t) zz E (t) 
+ Bap ee Hes u; (t) - v; (t) E 
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.(t) = u (t) [1 -u (t) - v (€) Im 


j-1 1,319 


+ [1 - u, (t) ιν. 
It should be noted that since the 
cov [X (€) , X4(t)] - E[X4 (€) X4(t) ] = E(X, (t) JE[X, (t) ] 


and since X,(t) is distributed N(0,0°); the Cov[X, (t) ,X,(t)] 
is equal to the mixed second moment, designated m, (t). 
Similarly, the variance of X; (t) is equal to the second moment 
of X,(t) - E[X,^(t)] = m, (t). 

Given q; (t) and m, (€) we have the distribution of Q; (t), 
since each Q; (t) is distributed Normally with a mean equal 
to 0,9, (t) and a variance equal to Qgm, (€). As a result of 
the possible correlation between ranks, of greater importance 
BONUS 2S the Multivariate Normal Distribution of Q(t). No 
routine was available that provided the capability of com- 
puting the probability of maintaining all the ranks within 
certain bounds; however, a bivariate normal distribution was 
used to compute probabilities of maintaining adjacent ranks 
within plus and minus one standard deviation from the expected 
rank level. The reason for desiring this capability will 


be discussed in Chapter VI. 
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V. BINOMIAL SIMULATION VS. THE DIFFUSION APPROXIMATION 


Prior to relying on any approximation of a real world 
system, a user should have a feel for how close the approxi- 
mation comes to reflecting the true system. In this chapter 
we will simulate the year to year operation of an imaginary 
Markovian 3-rank system. The system will be operating with 
Beemetion, attrition and accession rates such that the expected 
value of each rank is constant, i.e., the system is main- 
tainable at some specified level. The starting parameters 


are listed below, both in tabular and graphical form. 


- 


Q; (t) u; (t) v; (t) 

# of individuals Promotion Rate Attrition Rate 
Rank 1 1000 0.70 0.10 
Rank 2 2000 0.10 0.25 
Rank 3 1000 0.00 0.20 


The accession rate is specified as 0.20 of an available body 
of recruits estimated at 4000 individuals each year. 


Graphically the system would appear as depicted below: 


4000 





30 





The numbers in the parentheses are the expected number of 


individuals flowing across the appropriate arcs during each 


time period. Once again we note, that in a deterministic 


sense, the system 
leaving each rank 
the rank. 


Let 


2 


Q; (t) 


P. (t) 


ΠΣ 


u; (t) 


νι (5) 


a(t) 


is in exact balance, i.e., the number 


is exactly 


equal to the number flowing into 


total number in the system 


number of 


number of 
i at time 


number of 
i at time 


promotion 


attrition 


accession 


individuals in rank i at time t 


individuals promoted from rank 
t 


individuals attrited from rank 
E 


rate for rank i at time t 


tate for rank 1 at time t€ 


rate at time t. 


Note that P, (t) and W; (t) are not independent random variables 


since 


P,(t) + W,(t) < 0,(t) 


A program BIMOD, outlined in Appendix C, was designed to 


simulate this system. 


BIMOD operates as follows: 


E 





l. The program is designed to replicate the operation of 
a manpower system, for a user specified number of times. 
(Maximum number of replications allowed is 1000.) In the 
results which will follow, the number of independent repli- 
cations used was 300. A replication consists of 35 years of 
promotions, accessions and attritions of the manpower system 
under consideration. The time increment in each block, 
corresponds to "t" in all the variables mentioned above. 

The reason blocks of 35 years were chosen for comparison 

with the approximation, results from the fact that the 
diffusion approximation reaches a constant covariance matrix 
at the 35 year point. By terminating each replication at 35 
years, and storing the value of each rank at the end of each 
run, we can then compute an estimate for the covariance matrix 
of the simulation at the 35 year point. 

2. For each time period t, BIMOD generates a binomial 
random deviate with parameters ο, (Ε) and [u; (t) + v, (t)]. 
This deviate represents the total number of individuals 
leaving rank i at time t; both through promotion and through 
attrition. Call this random variable Li(t), A meets 

3. For each time period, we also generate independent, 
binomial random deviates with parameters L; (t) and 
U¡ (t)/10, (t) + v, (t)], ana O and a(t). These 
deviates represent the total number of promotions from rank 
i at time t, Pi(t); and the total number of accessions into 


the system at time t, A(t). 
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4. The Q; (t) are updated for each time period utilizing 


the following recursive equations 


Q, (t+1) Q(t) - L} (t) + A(t) 


Q; (t+1) Q; (t) - L, (t) + Pil (Ee fOr i = 2,3 


1 
In comparing the simulation against the diffusion approxi- 
mation, we would like to reduce the comparison to the random 
terms. We do this in order to determine how good the approxi- 
mation is in predicting the means, the variances, and the 
covariances of the random terms. 
As was discussed in Chapter III, equation (III.3) is the 


cornerstone of the diffusion approximation 


Dre O a "Qo Xi(t). (CDD 


We could describe the system of our simulation in similar 


terms, 


SE Q m. (t) 4 /Q. Ri(t) er) 


where Qm; (t) is the expected value of rank i at time t, 

M; (t) is the observed value of rank i at time t, and R. (t) 
is the random residual term for rank i at time t. The com- 
parison of (IV.1) and (IV.2) would commence with an examina- 


tion of the R; (t) residual terms 
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ο) -οπσπ, (5) 
R, (t) = O (IV.3) 
νο. 


O 


To see if the R; (t) do in fact exhibit an approximate normal 
distribution with mean equal to zero. 


BIMOD was run with the following input values: 


NUMBER OF REPLICATIONS 300 
YEARS CONSIDERED IN EACH REPLICATION 35 
INITIAL FORCE LEVELS 

Q1 1000 

Q^ 2000 

Q3 1000 
PROMOTION RATES 

u(t) 0.70 

u^ (t) 0.10 

u, (t) 0.00 
ACCESSION RATES 

a(t) 02 20 
ATTRITION RATES 

vi (t) 0.10 

v(t) 0.25 

va (t) 0-20 


At the end of each 35 year replication, the value of 
each rank was stored, and the residual values were computed 
as in equation (IV.3). The expected values used in equation 
(IV.3) were the initial force levels as stipulated above, 
since they are the levels we expect to maintain. When the 
300th replication was completed, the mean value of the residuals 


in each rank, ΕΒ.» was computed, and this in turn was used to 
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compute the estimate for the covariance matrix of the 35th 
year residuals. The formula used in computing the covariance 


matrix is 


(R. - R.) (R. - R.) 
Lp I ο 
Cov(Ri ,R4) D > 


where the R; are the stored values at the end of each repli- 
cation. 
The following is a compilation of the results obtained 


from the simulation. 


E[R; (t)] Kurtosis Skewness 
Rank 1 0.0398 2.7468 0.2054 
Rank 2 0.0162 2.6670 0.0248 
Rank 3 0.0041 2.8470 0.0708 


Kurtosis and skewness information was obtained from a 
histogram of the total run residuals. It should be noted 
that the Kurtosis Coefficient for a normally distributed 
random variable is equal to 3.0 and the skewness should 
equal 0. 

This first level of comparison of the R; (t) terms with 
the X. (t) terms seems to indicate that the assumptions of 
normality with mean equal to zero does in fact provide a 
reasonable approximation, with a possible exception of rank 


l, as indicated by its skewness coefficient. 
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The next level of comparison between R. (t) and X; (t) 
involves the predicted variances and covariances. Utilizing 
BIMOD the following covariance matrix was obtained for the 


R; (t) residual terms. 


Rank 1 Rank 2 Rank 3 
Rank 1 0552035 0.0000 2050022 
Rank 2 0.0000 0.4747 0.0245 
Rank 3 20.0022 0.0245 022733 


and this is compared to the following X. (t) covariance matrix 
obtained from the diffusion approximation program MNPWR, 


outlined in Appendix C. 


Rank 1 Rank 2 Rank 3 
Rank 1 0.2448 =0. 0210 -0.0005 
Rank 2 zo T 0.5059 "ο ο. 
Rank 3 -0.0005 20710289 0.2373 


At first glance, it would appear that the approximations 
for the covariance terms did fairly well. Of particular note, 
however, is the fact that the approximation specifies a 
negative correlation between ranks 2 and 3, while the simu- 
lation specifies a positive correlation for the same ranks. 

In order to better understand what is going on, the 


variance of the covariance terms, i.e., the variance of the 


E oen for all end of run Ri, was computed to 
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determine whether the diffusion approximation estimate falls 


in the range 9; Sl ο ο uelis The following 


are the variance results. 


Simulation 
Ranks Cov - 1 Std Dev Covariance Gov + 1 Std Dev 
2 -0.0173 0.0000 0.0173 
3 2201001562 -0.0022 0.0118 
23 0.0047 0.0245 0.0443 


While the covariance range specified above still implies 
a positive covariance for ranks 2 and 3, the covariance 
values under consideration are small enough to imply that 
the approximation estimate for the covariance is quite good. 

The next important matter of note, is the fact that in 
both matrices above, the covariance terms are essentially 
zero; and the fact that ranks appear to be independent of 
one another is borne out by very low correlation coefficients 
between ranks, both in the simulated system and in the diffu- 
sion approximation. 

A possible explanation for the resulting independence 
could be that since we do not vary promotion, attrition, or 
accession rates in either the simulation or the diffusion 
approximation, the probability laws of the models are the 
same as those in a model in which the motions of individuals 
through the system are independent. It is known that if 
the input random variables are independent and have Poisson 


distributions and the motions of the individual in the system 
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are independent, then the number of individuals in stage i 
at time t, σι (0), are asymptotically independent with marginal 
Poisson distributions as t > e [Refs. 1,8]. 

In comparing the variances, we note that the diffusion 
approximation provides reasonable estimates for the variance 
of the simulation residuals. 

In conclusion, one could possibly simulate a 6-rank sys- 
tem and obtain better results for the variances and covariances 
of the ranks. However, due to the enormous number of deviate 
generations required to produce representative output, it 
is felt that the much simpler computations of the diffusion 


approximation far outweigh their loss of accuracy. 
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VI. DIFFUSION APPROXIMATION RESULTS 


In this final chapter, we will present the results 
obtained when the 6-rank Supply Corps system is modeled 
utilizing the Diffusion Approximation Program, listed in 
Appendix C. In addition a program called MINI will be dis- 
cussed and its results reintroduced into a second run of 
the Diffusion Approximation to the Supply Corps system. 

The final subsection of this chapter will provide some 
recommendations for further study. 
A. A DIFFUSION APPROXIMATION APPLIED TO THE PRESENT LEVELS 

OF THE SUPPLY CORPS 

In Chapter II, we determined, at least in the determinis- 
tic sense, that under specific assumptions, desired rank 
levels, as defined by Ref. 9, could be maintained. We now 


wish to utilize those same attrition rates 


-. O02720. 1355,0.1121,0.0390,0.0790,0.1622] 


in the Diffusion Approximation and to apply them to present 


Supply Corps force levels. These levels are presented below: 


Rank Level 
ΟΞ 632 
Ο-2 691 
0-3 οι 
O-4 pP 
OE 616 
O-6 177 
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These values were obtained from the Supply Corps Officer 
Personnel section of the Naval Supply Systems Command, and 
were current values for 310X officers as of September 1979. 

The program MNPWR listed in Appendix C, incorporates the 
recursive equations derived in Chapter IV. A feature of 
this program is that, given the present levels, desired 
levels, and estimates for attrition rates specified pre- 
viously, the program will compute the following fixed promo- 


tion rates, 


E [0:5099,0.3711,0.1210,0.1005,0.0615,0.0000] 


and a fixed accession rate "a" of 0.0862. Note that these 
promotion and accession rates are those which are necessary 
to attain and maintain the desired force levels given that 
we start at the present force levels, and are the same as 
those specified in Section II.B of this thesis. 

Therefore, in the first of a series of computations, 
we attempt to see how long it will take to attain the desired 
force levels using the promotion and accession rates specified 
above. Note, we are not worried about maintainability, since 
Ref. l states for the deterministic case, that if a structure 
is maintainable, it is also attainable. However, in the 
stochastic sense maintainability has a slightly different 
connotation. In considering methods of control, we are now 
interested in methods that not only provide the expected 


rank levels, but also provide manpower planners estimates for 
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the variance of each rank. Returning to the initial compu- 
tation, the reason we are not concerned about maintainability 
is that if the program is able to compute non-negative 
accession and promotion rates, then the structure is main- 
tainable. 

A program run was made to determine how long it would 
take to attain the desired force structure. At a point 50 
years into the future, the expected number in each rank was 
approximately at the desired level. This length of time is, 
without much thought, unacceptable. A technique which allows 
faster convergence to the desired levels will be described 
in the next subsection. 

The following is a listing of the expected values and the 


covariance matrix of the system at time 50. 


EXPECTED VALUE OF EACH RANK AT TIME 50 


Rank E [Rank] Desired 
O-1 620.99927 621 
τ 624.99902 625 
055 994.99780 995 
Ο-4 862. 93999 863 
Ο-5 616.91138 617 
ους 233.85883 234 
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Analysis of the covariance matrix will be reserved until 


the next subsection. 


B. MINIMUM CONVERGENCE TIME 

The quickest method of attaining the desired force from 
the present force would be simply to promote (or attrit) 
all shortages (or excesses) in one time period. In the case 
being examined, this would be totally infeasible both for 
legislative and for practical purposes. To attain the 
desired levels in one time period would involve a very large 
number of officer promotions, which would be totally unaccepta- 
ble from the organizational standpoint. 

The approach taken in program MINI (see Appendix C), is a 
Dynamic Programming Approach. Given what are reasonable 
bounds for promotion and accession rates, the program will 
minimize the difference from present levels to desired levels, 
at each time period, for specified attrition rates. 


Utilizing the following bounds for accession and promotion 


rates, 
Lower Bounds Upper Bounds 

(accession and promotion (110% of the lower 
rates required to maintain bound rates) 
the system) 

acces- 0.0862 0.0948 

sion 

O-1 0.5099 0.5609 

ο-2 0-371 0.4082 

0-3 0.1210 0-1331 

0-4 0.1005 0.1106 

0-5 0.0615 0.0677 

Ο-6 0.0000 0.0000 
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and the attrition rates specified at the beginning of this 


chapter, we could reduce the convergence time required for 


the expected number of individuals in each rank to attain 


the desired level from about 50 years to at most 8 years. 


Rank 


O-4 
r5 
(9556 


As an output of program MINI, we obtain 


required to minimize the shortfall at each time period. 


CONVERGENCE TIME REQUIRED 


Time Required 
for Convergence 
years 
years 

year 

years 


years 


ο σι ο e U A 


years 


the promotion rates 


The 


following is a listing of the promotion rates suggested by 


MINI. 
TIME rl 
1 0.5099 
2 0.5099 
3 0.5104 
4 075286 
5 085292 
6 (082 92 
7 2 
8 025179 


MINI PROMOTION 
AND ACCESSION RATES 


0-27 

053902 
0.3754 
005829 
0.3904 
0.3904 
0.3904 
0.3904 


0.3790 


0=3 

0.1331 
021331 
0. I331 
0-1331 
39 3T 
0:1331 
0-1331 


0.1260 
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O-4 co ACCESSION RATE 
0.1106 0.0677 0.0862 
0: 11060.0677 0.0862 
027206020677 0.0862 
0.1106 0.0677 0.0891 
ΠΡ, 10) σσ, 0.0893 
0.1088 0.0677 0.0893 
0.1066 0.0677 0.0893 
ο Ἱοιο ὁ 0622 0.0875 





Note that between years 7 and 8 the promotion and 
accession rates begin to decrease. This is a direct result 
of the fact that we are approaching the desired levels, and 
following year 8, we then use u and "a" as specified in 
Section VI.A. 

Utilizing these computed accession and promotion rates 
as input parameters to the Diffusion Approximation yields 
information on how the system would react to system controls. 
The following is a summary of the Diffusion Approximation 


results at the 8 year point. 


Desired Expected Value 

Rank Level Approximation Result 
e-1 621 621.00903 

0-2 625 624.95386 

0-3 995 994.87744 

O-4 863 863.02441 

0-5 617 617.00854 

O-6 234 233.97964 


Of even greater significance are the variance/covariance 
results. Below is the covariance matrix obtained at the 8 
year point. Variances for each rank are on the main diagonal. 

A striking result becomes evident if one compares the 
expected value and the variance of each rank. The indepen- 
dent motions property is still evident, even with promotion 
control, i.e., the expected value in each rank is approximately 
equal to the variance of the rank (as one would expect for 


Poisson distributed random variables). What this may imply, 
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is that in this case even though we vary the promotion and 
accession rates, the changes allowed in those rates are 
much too small to counteract the independent motions property. 
As a result, bivariate probabilities taken indicated that 

the ranks are essentially independent. 

If in fact the ranks are relatively independent, one 
could then compute expected value plus and minus one standard 
deviation of force level in each rank, and know that with the 
assumption of a normal distribution for force levels, we 
have approximately a 68% chance of being in that interval. 
From our calculations we get the following results for upper 


and lower bounds of force level. 





Rank Lower Desired Upper 
Bound Level Bound 
SA 596 Sex 645 
9-2 598 625 651 
03 962 995 1027 
0-4 833 863 892 
0-5 SoL 617 642 
0-6 218 234 249 


Given these values, manpower planners need not overreact 
to random deviations in force levels within the above bounds. 
If, for example, the O-4 rank were to fall to as low as 833, 
this would not necessarily imply that we should increase 
next year's accessions; it could be rationalized as a reason- 
able value as compared to the variance computed by the Diffu- 


sion Approximation. Without the variance estimates, planners 
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do not possess the tools necessary to determine whether or 
not fluctuations in each rank can be explained by the internal 
randomness of the System, or whether the fluctuations are 
due to external forces which require changes in accession 


and/or promotion rates. 


C. RECOMMENDATIONS FOR FURTHER STUDY 

The Diffusion Approximation hinges on reasonable estimates 
for the attrition rate of each rank. Therefore, Recommenda- 
tion 1 would be to derive a useful estimating relationship 
for the forecasting of attrition rates. 

Another possibility for an extension of this thesis, 
would be to consider confidence intervals for attrition 
rates, and to input these intervals with the Diffusion Approxi- 
mation. Analysis of the results obtained would be interesting 
in regards to the independent motions observed when point 
estimates were used for attrition rate. 

It is clear from the foregoing analysis that maintaining 
a grade structure in a stochastic environment is much less 
straightforward than a deterministic analysis of the problem 
would suggest. Under rather special conditions we have seen 
that it is possible to maintain a structure in spite of random 
variation in each rank. With better forecasting of attrition 
rates, it is felt that the Diffusion Approximation will prove 


to be an invaluable tool to manpower planners. 
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APPENDIX A 


DESIGNATOR 3100 
AUTHORIZATIONS SUPPLY CORPS (SC) 


IVA 


ει 
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APPENDIX B 


DERIVATIONS 


B.1 DERIVATION OF THE q.(t+l) RECURSION, THE VARIANCE 
RECURSION, AND THE RELATIONSHIPS FOR Se 


Q,(t+1) -Q,(t) -» a(t)Q, + Vatt) [1-a(6)O w ^ (t) 
- u, (t)Q, (t) - Ya, (ey [1 - a, (10, (€) w,^ (t) 


- v, (t)Q,(t) - “vp Ce) Tv, (6) 10, (ey νι" (5) 


(for definition of the mentioned variables, see Chapter IV). 
Substituting YQ. xX, (t) + Q.q,(t) for Q,(t) and YQ, X, (t+1) 


T 9,4, (t+1) FOr οι (t+1), we have 
BN EH) -X,() * QS ay (t+1) - 40, ay (t)1 = 


07140, att) + vaii alt Wo tt) 


O 


EU EX (t) “ο αι (Ε)] 


- afu, (t) [1 - 0, (6) 101//8)X, (€) + a) (t) wW” (t) 


- vQ (C) D, (€) + YQ. q, (t)] 


-Av, CO. v, CO IGLV9) X, (€). 9 a1 (01 W^ (0) 
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X, (t+1) > X, (t) = 


Ὡς 9, (6) - Qf a, (ttl) + VO alt) + Vatt) [1-a(€)] w^ (t) 


uj(t) Xi(t) - ui(t) YQ, q, (t) 


- Ju, (€) [2 7 u (E) ILQ/40) X, (€) & q, (1. wj? (6) 


E, (EY * v.(t) YO. q. (t) 


Av, (t) (1 - v4 COTEGA/QO X (t) +a (1. wy Ct) 


Heuristically, since Xj (t) is a normal random variable, 


X4 (t1) - Xi (t) must be finite and therefore any terms with 


νου must approach zero as Q > œ. Therefore in the limit as 


> © 
Qo 


X, (t+1) - X, (t) πι ach WP (t) 


A Ἑ 
uy (4) [1 uy (t)1q, (t) W (t) 


E u(t) X, (t) 


E MACULIS CEN ee i 
y, (€) X, (€) - Vv, CO [1- v, (€)1q, (€) W,^ (€) 


Rewriting, we have 
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X, (t+1) = [1 - u, (t) - v, (t)] X, (t) 


+ vai- at] w’ (t) 


πο Ta τα E 
u (t)i-uj(t)aq(t) WW (t) 


——  — — ας E 
"v, (€) [1 2 v4 CE) Ja, (€) Wi (t) (B.1.1) 


It also follows that in the limit as Qo > 0 
q, (t+1) = qi (t) = a(t) - tu, (t) + νι (8) αι (5) 
or 
ΠΕΙ - αἰεὶ: [1 - u, (t) - v, (t)]Ig, (t) UBL) 


Using a similar approach for Q, (ttl) - Q; (t), i # l, the 


following result was obtained in the limit as Q, > o, 
X; (t*l) = αι (8) νι (8) + [1 - u; (t) - v, (t)] X; (t) 
πηι [ΕΠ (€) π. ‚Fo 
τ "τι icd 
/.(t)ii-u.(t)]aq. (€) W,^ (t) 
- BM, ΥΠ i 


οτι gs . 9 9 L . ο 
V7, (ey tL - v, (6) Ta, (E) Ww," (t) 3) 
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and 
q; (t+1) - q, (t) = u; (©) qd; (t) = [u; (t) +v, (t) Iq, (t) 
ΟΥ 


q; (t+1) = u_, (tq, (8) + [1 -u, (t) - v, (t)]q, (t) 


(B.1.4) 


We now have expressions for all the X; (t*1) and αι (511) 
in terms of X; (t) and q, (t) respectively. 

In order to calculate the variances and covariances 
of the {X,(t),i=1,...,r}, we must first determine the first 


and second moments of the X, (ttl). 


E[X, (t+1) ] = [1 - u, (t) - v, (t)] E[X, (t)] 


+ va(t){l-at)) E(W.*(t)] 


Va, (ey (1 a, Ce) 1a, Ct) ElW,* (t)] 


UT OT (0 sm” (t)] (B.1.5) 


Recall that in oe and λα ee eee 
are normally distributed with a mean of zero. The random 
variable X4 (t) represents the randomness in the system at 


time.t. Initially, since the system is starting from a known 


S 





set of rank levels, X4 (0) ls defined to be zero. Asa 


BouPt, equation (B.1.5) implies that E [X] (t+1)] πο τος -all 


time t. Further, 
Var[X, (t+1)] = [l-u (t) sg Reg Var [X, (t) ] 
ΠΠ ae Var[w ^ (t)] 
+ 9, (t) [1-u, (t)1q, (€) Var(W,* (€)] 


+ v, (t) [1 - v, (6) ]q, (6) Var [W," (t) ] 


Since the variances of the Wi random variables are equal 


to one, 
Var[X, (t*1)] = [1-u, (€) - v, (€) 1^ Var [X (t) ] 
+ a(t) [L-a(t)] + ui (t) [1 - uy (t) Iq, (t) 


+ v, (t) [1 7 vi (€) la, (€) 


Since the expected value of X, (t) is equal to zero, the 
variance of X, (t) is equal to its second moment which we 


will call m, (t). Therefore we have the following: 


m (tel) = [l-u (t) -v, (01° m, (t) + a(t) (1-a(t)] 


+ tu, (t) [1 - u, (4) ] + v (t) [1- v, (t) ]9 ai (€) 
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Using similar reasoning, we get a general recursive scheme 


FOr 1i > |: 


m, (t+1) = πο πο 


TET 
"EINEN E N E) 

ak 1 d 
d EIE ae 


* tu, (t) [1 - u, (t)] + v¿(t) [1 - v,(t)])q, (t) 


a t) 


+ 21u,_, (8) [1 - u, (t) - v, (8)]) Di-l,i 


Note, however, for i» l we get a covariance term. In 
the next section, we will obtain a recursive scheme for 


computing the covariances. 
B.2  DERIVATION OF COVARIANCE FORMULAS 
Recall: 
X, (t1) - [1 - u, (t) - v, (t)] X, (t) 
E 
+ va(t) [17a(t)] W, (t) 
- a Œ) ll- u, (t) jq, (©) W” (t) 
l l l l 


ED a Ea mb 
v, (t) {1 vı (t) lq] (t) W ΓΕ, 


55 





If we let 


bi (t) Ξ 1 - u; (t) = v; (t) 


ae) ie att)! 


ee) 


h, (t) = u; (t) [1 - ui (t)] 
BL INDE EI 
then 
Xj(EM) - b,(t) X,(t) + vote) Wo (t) 


πες το Wee επ... 
hi(t)qi (t) W (t) k4 (t) q4 (t) Wi (t) 
Similarly recall 


Χρ (ΕΓ) = u (t) X, (t) ~ b, (t) X, (t) 
E oa o e a = vha(t)q.(t) WT (t) 
l T l 2 2 2 


πας (EY Ww o 
SEL WA (t) 
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g 


yy 


>) ¡Y " (3) ,U -1 * LUI 


= 
o 


— 1356 - Diss = 4312 


"x 


It is assumed that at time zero the Cov[X, (t),X,(t) ] = 0 


feet = O0, 
Cov[X, (t+1) ,X, (ttl) ] = E[X, (t+1)X,(t+1)] - ElX, (t+1) ]E( xX, (t+1)] 
x (t+1)%, (ttl) = 
[b, (t)Xx, (t) + Yc(t) w^ (t) — vh, (t)a, (€) Wt (t) - vk, (t)q, (t) WU (t)] 


᾿ Ρ E P 
[u, (€) X, (€) +b, (t)X, (t) + vh, (t)q (t) W (E) vh, (t) q, (t) W, (t) 


- AO IO Ww, (£)] 


2 2 P 

= bi (t)u (€) X (t) +b, (t)b, (€) X (t) X, (t) * bi (t) /h (t) q, (€) X, (€)W, (t) 
P L 

- b. (t) ^h. (£) a, (£) X, (E) W, (t) - b- (t) vk, (t) a, (E) X, (€) W^ (t) 


P P 
+ u, (t) vc(t) X (tW (t) + b, (t) va (t) X2 (t)W. (t) 


+ 


Tape Ty wP EN = an P 
(Eh, (6), (E) 1, (6), (t) - ve(t)h(E)q, (E) W, (tW (t) 


ae we P 
c(t)k, (t) 9, (5) W, (t) W (t) 


u, (t) “Ay (ta, (6) X, (€) W," (t) - b,(t)/h Tg, le) X, (€)W," (€) 
- h, (t)q, (t) [W Pity}? + vh, (t)q. (t)h, (t) q, (t) W P(e)w Pt) 
d; 1 Ego l 
ann PME L 
* /h; (Eq, (E) k, (E), (E), "(E)W,^ (€) — 0, (6) vk (ea EX} (tw) (8 


Τι Ῥ E 
- bj(t)/k (jq; (€) X, (00,78) - qj (€) Vh; (EK, (€) wj" (tw (t) 
+/K Tea Th UOS NW. dw Pee) ka tkigse 9 wPew, P 
1 (699, (8)h, (4) 9, (EW, 1 1 (6) q, (t) k, (t) qp 2 1 
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ElX, (t+1)X, (t+1) ] = 


2 p 
b, (€) u, (E)E[X, “(t)] +b, (t)b, (t)E[X] (t)X, (t) ] +b (t) VA, (E) q, TOELX, (630,5 (6) 
- b, (€) 'h; (Oa; (OEIX, UA (t) ] - b, (€ "E; (Oa (B EX, (tw, (t) ] 


+ u (t) Ve (EJEIX, (€)W. (t)] * b, (t) /c(E)E[X, (yw P (t) 


wc(E)hy Eg) EI, ()w ^ (€)) -Vc(OR, Og, (0 ED (t)W_ (t)] 
τετ n P E P 
c(E)k,(E)q, (E) ED, (e)w ^ (€)] - u, (6) Ya, TEYq, (E) Elx, (6), (8) 


- b, (t) ^h. (t)g, (t) E[X, (t)W," (t)] = h, (t)q, (EEL (H,” (t)) 73 


+ 


Yh, (e)q, (ein, Wa, (e) EM, (6), (t)] 


L L 
u, (t) Vk, (t)q, (E)EIX, (t) W) (t)] - b, (t) (€) a, (E) E[X), (t) W) (E) 


P L Παπ; E L 
αι (5) h. (€) Κη (€) E [W] (t) W ΓΕ 1 3, (t)h, (t)q, (t) ElW, η (EJ) 


+ 


/k (Oa (Ok; (Og; Elm) (e), (t) ] 


Since the X; (t) are assumed independent of the W; (t) 


and EL" (6) *] = 1, the above simplifies to 


m, 4(t*1) Ξ E[X, (t+1)X, (t+1) ] = b, (t) uy (t) m, (t) 


+ b, (t) b. (t) m,..(t) - h4 (t) q, (t) 


1994 
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my 7 (t+1) = [1 - u, (t) - V(t)] u (t) m, (t) 


κ | -u (t) -v (t)] ils u, (t) E v4(t)] my 2 (€) 


- 0, (t) [1 -u (t)] q, (t) 
A similar approach can be taken for all cross product 


terms, and the following general relationships will be 


obtained: 


me (t+1) = ΘΙ [1 - u (t) -vQ (9l πο 
+ [1 - u, (t) gr ee ur m, (8) 
NET, )- itl 
Det) = u; Q(t) u; (t) mi 1,1 € 
+ u; Q(t) (1 RED mr a 


+ u; (t) D! -u,(t) λε m, (t) 


* [1-7 u, (t) E 7 ο, ie ae) 


u, (t) [1 - u; (t)] qj (t) 
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> l, j £ i+l 


u = u, ,(t) Den ΙΙ 


ας (τα) τν (ο πι ¿(t) 
+ e [1 - u, (t) - v, (t)] m je 


+ {l-u, (t) -v, (t)] ο a πο e 


To fill out the remainder of the covariance terms, recall 


that 


t E Mes VELA : 
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(2)+0(3) 


RUN 
CCEIR(DSEEC,CZERO, A1) 
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BISCMIAL SIMULATICN PROGRAM (FROGRÁAM EIMC[) 
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